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Abstract 

Bilinear forms for some nonlinear partial difference equations(discrete soliton equations) are de- 



' rived based on the results of singularity confinement. Using the bilinear forms, the Af-soliton and 



algebraic solutions of the discrete potential mKdV equation are constructed. 

Keywords: discrete soliton equation; bilinear form; Hirota's method ; singularity confinement; dis- 
crete potential mKdV equation 

1 Introduction 

Hirota's method is one of the most powerful methods to study soliton equations. The crucial step in 
the method is to find the proper dependent variable transformation (DVT) which reduces a given soliton 
equation to bilinear (or multilinear) forms, that is, to express the dependent variable in terms of a t 
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function. Once we get the bilinear form, we are able to construct systematically particular solutions 
including the multisoliton solutions. 

The bilinear approach is also useful for discrete soliton equations. In order to obtain the integrable 
discrete versions of soliton equations, Hirota discretized the bilinear forms of the soliton equations so as 
to preserve solution structure^ |). Then finding a proper DVT, we are able to construct the nonlinear 
difference equation. 

But when we try to find the solutions for a given nonlinear difference equation, it is difficult to find a 
DVT which reduces a given discrete soliton equation to bilinear (or multilinear) forms, even if we know 
the results of corresponding continuous case. 

Recently, Ramani, Grammaticos, and Satsuma proposed to use the singularity confinement test to 
get proper DVT's for nonlinear ordinary difference equations, and derived those for the discrete Painlcvc 
equations ||. Singularity confinement (SC) was introduced as the discrete analogue of the Painleve 
property for continuous systems, that is, as a criterion for integrability of discrete systemsQ. In this 
article we apply this method based on SC to some nonlinear partial difference equations to find the proper 
DVT's and the bilinear forms for them. In particular, for the discrete mKdV equation we construct the 
./V-soliton and algebraic solutions. 

2 Discrete KdV Equation 

We first consider the discrete KdV equation of the form^) 

(4 -4)- (l> 

We have a singularity whenever the it™ in the denominator happens to vanish. This has as a consequence 
that both u™ +1 and u™_jf diverge, whereupon vanishes again and u™^ 3 , u™_j~ 2 3 ; u ™-2 are finite. Thus 
the singularity is confined. Following the SC criterion, therefore, the discrete KdV equation is integrable. 
Singularity pattern of this case is 

{C,< + \<-^<i 2 } - {0,co,co,0}. (2) 
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We use this information in order to express u in terms of r functions. Let us assume that the singularity 
pattern is caused by the fact that a r function t™ becomes zero at (m, n). Then, the singularity pattern 
suggests the transformation 

mm— 2 

< = ~ i ■ (3) 

n m—1 m—1 v ' 

~n r n+l 

Since and u™.^ 1 contain in the denominator, and u m and contain in the numerator, 

eq.(|^) is certainly consistent with the singularity pattern. Substituting eq.(|^) into eq.(Q) we obtain readily 
the trilinear form 

r m+1 m— 2 m— 2 m— 1 m m— 2 c m— 3 m m m — 1 m m — 2 /,|\ 

ot ti T n+1 T n+2 er ri T n+l T ri+2 — 0T n+2 T n T n+1 tT n+2 T n T n+1 • V*J 

Muliplying eq.([|) by t™^ 1 and manipulating the result, we obtain 

£_m+l m-2 _ m-1 m 
nm _ pm-1 pm _ "'ft 'n+1 c 'rt /r\ 
r n — r n+l ) e n — m-1 m ' ^ ' 

T n+1 T n 

This gives rise to F™ = a{m + n), here a(m + n) is an arbitrary function of m + n. Thus we obtain the 
bilinear form: 

^C +1 C+1 2 - eC-V- +1 - a(m + njr^r™ = 0. (6) 

In case of a(m + n) = 1, this bilinear form is the same as that obtained by Hirota|^|. From this we can 
obtain the solutions through formal perturbation procedure. 

3 Discrete To da Equation 

Secondly, we consider the discrete Toda equation of the form[|) 

rm+1 rrn T/m T/m+1 (n\ 

1 n 1 n — v n v n-l ' \' J 

fm+l|/m+l rm T/m /o\ 
L n V n ~ 1 n+l V n • (°) 

We have a singularity whenever the V^™ happens to vanish. The singularity pattern is 

{VZW^V^ 1 ^ 1 } - {0,oo,oo,0}, (9) 
{C+i.C+iSC^C^ 1 } - {co, 0,0, co}, (10) 
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and the singularity is confined. Let us assume again that the singularity pattern arises from that a r 



function r™ 1 vanishes at (m,n). Then, the expressions for V and J, dictated by the singularity pattern, 



are 



m m+1 

T n~l T n-l 
m+1 m 

jrn = 'n-l n-2 (12) 
T n-l T n-2 

Substituting these into eqs.(0) and (||), we obtain readily 

1 - C+iVr-V = «M[C] 2 - (13) 

Here, a(m) is an arbitrary function of m. In case of a(m) = 1, the above bilinear form is the same as 
that obtained in Ref.jjJ. 



4 Discrete Hungry Lotka-Volterra Equation 



Next, we consider the discrete Hungry Lotka— Volterra equation of the formj 



M -, , r t 

TT 1 + (u) 



n i=l 1 ""n+i 

Examining eq.([i"4|) in the form solved with respect to u^^ M , we find the singularity pattern 

{ U l-M-1> "!,+!/} {0,CXD, OO,0}, (15) 

and the singularity is confined. However, when we assume the singularity pattern arises from vanishing of 
a t function, we can't obtain the expression of u in terms of the r function which satisfies the singularity 
pattern. So we divide the singularity pattern into two patterns: 

{«n-M-U u n-l' u w «n+M> ~> {0. °0, 0}, (16) 
{CiX-n^ 1 ."!.} -» {oo,-i -i,oo}. (17) 

These patterns suggest the transformations 

r t+i _t 

t _ 'n+M+l'n-M ( , R \ 

a n — _t+l t \ 10 ) 



'n+1 'n 

-t _t+l 
'n+1 <n 



T n+l T n 



-, (a : const). (19) 
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We thus obtain the bilinear equation of the form 

KIm+A-m ~ a§< +1 ri +1 + rl+Vn = 0, (20) 
which is the same as that obtained in Ref. 0] . 

5 Discrete Potential mKdV Equation 

Let us next consider the discrete potential mKdV equation of the form|| 

nv m+1 - v m+1 

v rn+2 = ^ -1 ^ (21) 



To apply SC, it is convenient to consider 



jm+l _ Tm-l r 1 n+l [ ll n 1 
h L1 n+l 1 M 1 n 



(22) 



where 7™ and are related as I™ = v™/v™_i. In this case we have two singularity patterns. Either 
I™ first takes the value l//x, in which case we have 

{C, C + \ C-Y. C- + i 2 } ^ 0, oo, fx], (23) 

A* 

and the singularity is confined, or we start with /i, in which case we have 

{C.C+SC-V-C+i 2 } -> {/wo i}, (24) 



and the singularity is confined. We assume that the singularity patterns (|23| ) and ( [24|) come from vanishing 
of r functions and G™ at (n, m + 1). Then it is suggested that J™ are in the forms 

pm-l 
J™ = n + -H±J-P 
r n+l 
1 prn+1 

= - + — Q 

" R, (25) 



rprn 



-i /^m — 1 
_ 1 , U n+1 rp 



1 1 

-U, (26) 
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where P, Q, R must be expressed in terms of G™'s and S, T, U in terms of F™'s. Combining eqs.(|25|) and 
(p6|), we find the following simple expressions for I™ 

4 = V + a pm rm (27) 

-i jpm+lrrm— 1 

= " "+ 1 f29) 

where a, (3 are arbitrary constants. The bilinear equations are obtained from the equivalence of the above 
three expressions of J™ . 

,.pm/mi Trim /OWl a nm+1 sym— 1 ("qi^ 



The equation fl22| ) is, then, satisfied automatically. From the expression fl29|), we have 

m = U n+1 

n fP m 
r n+l 

In order for this to satisfy (plf), f3 = —a//i is required. 



(32) 



6 Discrete mKdV Equation 

We finally consider the discrete mKdV equation proposed by Tsujimoto and HirotaQ], 

w l k +1 (l + 5w l k +\) = + (33) 
1 + aw 1 ^ 1 l + aw{ ' 

Bilinerization of eq. (|33|) is done in a similar way to the previous section. We have two singularity patterns, 
either w u takes the value — - in which case we have 

{w l k ,w l k+1 ,w l k %w l + + \} -> {-i -~ oo,0}, (34) 

and the singularity is confined, or we start with 0, in which case we have 

M V^W+a } -» {°> «>> - J > ™ }> (35) 



G 



and the singularity is confined. From these patterns, we obtain the expression for w k , 

w l k = + a — — , (36) 

k 1 - 1 a 1 

= P J=i a l ■ ( 3? ) 

K fc-l°fc-l 

1 K Z rr i-1 

= ~7 +^ i-1 I ' 38 ) 

where a' , (3' and 7' are arbitrary constants. Then we get the bilinear equations from the equivalence of 
these expressions, 

5 7 '4 +1 4 - 44 +1 - /SWfc V+i = , (39) 
( 3\ J „i+i 'slJ+^J 



H*h °k -\^-- a ) k k a'KTi^-i = . (40) 
Comparing these equations with eqs.(|30|) and (pl|), then we easily find that they are equivalent if we put 

m = —k, n = l, F = «, G = <r, (41) 

and 

Finally, we obtain the explicit relationship between discrete potential mKdV equation and discrete mKdV 
equations. Noticing the relation eqs.([ll]) and (|42|), we find from eq.(|3|)that w and v are related as 

--K^f- 1 )' (43) 

which can be also verified directly. The point is that it requires "art of discovery" to find the above 
relationship just by looking at the equation, but now it becomes a "routine work" if we look at the 
bilinear equations which is obtained through SC. 



7 A/-soliton Solution of the discrete potential mKdV Equation 

Let us construct the N— soliton solution of the discrete potential mKdV equation. Indeed, we can use the 
perturbative technique to obtain the iV-soliton solution, but here we derive it from the reduction from 
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the discrete KP equation. The bilinear form of discrete KP equation is written as 



a 1 (a 2 - a 3 )r(ni + 1, n 2 , nz)r{n\, n 2 + l,n 3 + 1) 
+a 2 (a 3 - ai)r(ni,n 2 + l,n 3 )r(ni + l,n 2 ,n 3 + 1) 
+a 3 (ai - a 2 )T(n 1 ,n 2 ,n 3 + l)r(ni + l,n 2 + l,n 3 ) = 0, 



(44) 



where r depends on three discrete independent variables n^, n 2 and n 3 , and ai,a 2 and a 3 are the difference 
intervals for ni,n 2 and n 3 , respectively |J. The Casorati determinant solution to eq.([44|) is given by |TT] 

fi(n 1 ,n 2 ,n 3 ;s) /i(ni, n 2 , n 3 ; s + 1) •■• /i(ni, n 2 , n 3 ; s + N —1) 
/ 2 (ni, n 2 ,n 3 ;s) / 2 (ni, n 2 , n 3 ; s + 1) ••• / 2 (n x , n 2 , n 3 ; s + JV - 1) 



r(m,n 2 ,n 3 ) 



/iv(ni,n 2 ,n 3 ;s) /jv(ni, n 2 , n 3 ; s + 1) ••• n 2 , n 3 ; s + iV - 1) 

where fj's are arbitrary functions of ni,n 2 and rt 3 which satisfy the dispersion relations 

^n k fj(ni,n 2 ,n 3 ; s) = /j(rii, n 2 , n 3 ; s + 1) (fc = 1,2,3). 

Here A nk are the backward difference operators: 

A nk f(n k ) = (k = 1,2,3). 

For example, we can take /j as 

/iCni.na.najs) = pj(l - ^ai)"" 1 (1 - P^)"" 2 (1 - P^s)"" 3 

+q s j {\ - 9j ai)-' ll (l - %a 2 )-" 2 (l - ^a,)-" 3 , 

where j>j, qj are arbitrary constants, then it will give the A^-soliton solution. 

Now we assume that there exist a non-zero constant $ such that for arbitrary ni,n 2 , n 3 [ft0| 

r(ni,n 2 ,n 3 - 2) = $r(ni,n 2 ,n 3 ). 

Then defining as F(ni,n 2 ) — r(ni, n 2 , n 3 ), G(ni,n 2 ) — r(ni,n 2 ,n 3 + 1) we have 

a 1 (a 2 — a 3 )F(ni + 1, n 2 )G(ni,n 2 + 1) 
+a 2 (a 3 - ai)F(ni,n 2 + l)G(ni + 1, n 2 ) 
+a 3 (ai - a 2 )i ;l («i + l,n 2 + l)G(ni,n 2 ) = 0, 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



and by shifting n 3 — > n 3 + 1 in eq. (p4|) and using eq.(fl9| 



01(02 - a 3 )G(ni + l,n 2 ).F(n,x,n 2 + 1) 
+a 2 (a 3 - ai)G(m,7i2 + l)F(n\ + l,n 2 ) 
+03(01 — a 2 )G(ni + 1, n 2 + l)F(ni,n 2 ) = 0. 



(51) 



After shifting n% — > nx — l,eqs.(50) and ( pl| ) are transformed into the following equation 



a x (a 2 - a 3 )F™G™ +1 + a 2 (a 3 - ax)F™ X G™ = a 3 (a 2 - a 1 )F^ 1 1 G^ +1 , 
a 2 («3 ~ a{)F™G™ +1 + 0l (a 2 - « 3 )C + iC = « 3 («2 - ai)iC +1 G™+x\ 



(52) 
(53) 



through the variable transformation 



m 



-m - n 2 , 



n 2 , 



(54) 



where F(ni,n 2 ), G(ni,n 2 ) arc defined by 



F(ni,n 2 ) = F? = F-^- n \G(n u n 2 ) = G™ = G~ 2 ni " 



(55) 



Moreover if we must choose as 

_ a 2 (a3 - ax) _ a 3 (a 2 - «x) „ _ «3(«2 - «x) 

ai(a 2 -a 3 ) ai(a 2 -a 3 ) a 2 (a 3 - ax) 

eqs.(|30|) and ([n]) reduce to eqs.([52|) and (|5 



(56) 



Now we impose some constraint on the parameters of the solution so that the condition (49) is satisfied. 
We observe that 



fj(ni,n 2 ,n 3 - 2; s) 



= P °{1 - Pj a x )- n -{l - Pj a 2 )- n -{l ~ P] a 3 )- n ^ + gj(l - <&ai)- ni (l - q,a 2 )-" 2 (l - q 3 a 3 ) 



\-n 3 +2 



= p s j(l - Pjai )- ni (l - Pj a 2 )- n *(l - Pj a 3 ) 



-n 3 +2 



1 - qja 2 \ I 1 - Oja 3 \ / 1 - gjaa 



1 - Pj a 2 



1 -pj-a 3 



1 - Pj a 3 



(57) 



Thus, if we put 



1 - a 3 pj = -(1 - a 3 qj) , 



(58) 



Qj = — -Po 

^3 



(59) 



then wc have 



f j {n 1 ,n 2 ,n 3 ~ 2; s) = (1 - Pja 3 ) 2 fj{n 1 ,n 2 ,n 3 ; s) 



(60) 



Using this relation we readily find the following relation: 

JV 



T(n 1 ,n 2 ,n 3 - 2) = ]^[(1 - p k a 3 ) 2 T(n 1 ,n 2 ,n 3 ) , 



(61) 



fc=i 



that is, the condition ( |49| ) is satisfied for this solution. 

Thus through the variable transformation ( fii] ) and (fio]), the iV-soliton solution of discrete potential 
mKdV equation is given by 

... 



r n+l 



T\ 



n.+ l 



4>i(m, n;s) <j>\ (to, n; s + 1) 
4> 2 (m,n]s) 4> 2 {m,n;s + 1) 



</>i (m, n; s + TV — 1) 
(/> 2 (m, n; s + iV — 1) 



^Ar(m,n;s) </>jv(m, n; s + 1) ••■ 0at(to, n; s + TV — 1) 



G7", 



ipi(m,n;s) %l)i(m,n] s + 1) 
ifj 2 (m,n;s) ip 2 (m,n; s + 1) 



t/)i(m, n; s + N — 1) 
ip 2 (m, n;s + N—l) 



^at(to, n; s) ^w( m : n ! s + 1) '"' ^jv(to, n; s + iV — 1) 



^■(m,n; s) = ^(1 - Pj ai) m+n (l - pjOa)"" + #(1 - <Z,ai) m+ "(l - fc-aa)" 
^•(ni,n 2 ;s) = p. s (l - p j0 i) m+n (l - ft-aa)"" - - <?,ai) m+ "(l - ^a 2 )" 



9j 



«3 



For example, the 1-soliton solution is given by 

y7n = Pl(l - P iai) m+n+1 (l - Piaa)-"- 1 - gf(l - g iai) m+ " +1 (l - gi^)-"' 1 



(62) 



(63) 



(64) 



(65) 
(66) 
(67) 



(68) 
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8 Algebraic Solution of the Potential mKdV Equation 

Let us construct the algebraic solution of the discrete potential mKdV equation through the reduction 
of the discrete KP equation. 
First notice that 

d 

f j (m,ri2,n3;s) = — p^l-pj-oi) ni {l-p 3 a 2 ) n2 {l-p 3 a 3 ) ns , 

Opj 

= -+Et3t ^-¥ir(i-ft«ar(i-ft»sr > m 

\ p i k=i p i ak J 

satisfies the dispersion relation (Effl) , and hence this also give the solution of the discrete KP equation 
(ff4|). Moreover, if we notice that the discrete KP equation (Q) is invariant under multiplication of 
c iQi 1 C 2Q 2 2 c 3 % 3 on r function, where Cj and q 3 , j = 1,2,3 are arbitrary constants, we find that it gives 
the algebraic solution of the discrete KP equation. More generally, fj can be chosen as 

f 3 im,n 2 ,n 3 ) - ^p^l - p jai )- n '(l - Pj a 2 )- n2 (l - Pj a 3 )- n3 , (70) 
dp/ 

where kj is an arbitrary integer. In particular, we define the polynomials by 

oo 

E^(ni,n 2 ,n 3 )A fc = (1- Aa 1 )-" 1 (l-Aa 2 )-" 2 (l-Aa 3 )- n3 (P k = 0, for k < 0). (71) 

fc=0 

The simplest examples of Pk are given by 
Po = 1, 

Pi = <2i7ii + a 2 n 2 + a 3 n 3 , 

P 2 = ^{ai?ii + a 2 n 2 + a 3 n 3 ) 2 + ^(afni+aln 2 + aln 3 ), 



Then, the r function 



T N (ni,n 2 ,n 3 ) 



P n (ni,n 2 ,n 3 ) 
Pi 2 -i(n 1 ,n 2 ,n 3 ) 



P n+1 (ni,n 2 ,n 3 ) 
P i2 (m,n 2 ,n 3 ) 



Pi N -N+i(ni,n 2 , n 3 ) P iN _ N+2 (ni,n 2 , n 3 ) 



Ph+N-x{ni,n 2 ,n 3 ) 
Pi 2 +N-2(ni,n 2 ,n 3 ) 

P iN (n 1 ,n 2 ,n 3 ) 



(72) 
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where i± > ii > ■ ■ -ijv > are integers, is nothing but the Schur function solution of eq. ( f44j ) labeled by 
the Young diagram i 2 , ■ ■ ■ , «jv)- 

In order to satisfy the condition (E9J) , let us consider the r function, 



T N (ni,n 2 ,n 3 ) 



P N (ni,n 2 ,n 3 ) P N+1 (n 1 ,n 2 ,n 3 ) 
PN-2{ni,n 2 ,n 3 ) P N -i(ni,n 2 ,n 3 ) 

P-N+2(ni,n 2 ,n 3 ) P^ N+3 {n\,ri2,n 3 ) 



P2N-i{ni,n 2 ,n 3 ) 
P2N-z{n\,n 2 ,n 3 ) 

Pi(ni,n 2 ,n 3 ) 



and impose a condition, 



Pk(ni,n 2 ,n 3 - 2) = P fc (ni, n 2 , n 3 ) - c 2 P k - 2 (ni,n 2 ,n 3 ), 

where c is some constant. Then it is easily verified that the r function satisfy the condition (ff9|). 
For simplicity, we introduce a function, 

9( ni ,n 2 ,n 3 ) = (1 - A ai )- ni (l - Aa 2 )-" 2 (1 - Aa 3 )~" 3 , 

which is a generating function of the polynomials Pk- Then eq.(|74|) implies that 

2 \2, 



(ni,n 2 ,ra 3 - 2) = 0(ni,n 2 ,n 3 ) - c A 6(ni,n 2 ,n 3 ) , 



from which we have 



1 - Aa 3 = ±(1 - c 2 A 2 )5 



(73) 



(74) 



(75) 



(76) 



(77) 



Defining as F(ni,n 2 ) — T^(ni,n 2 ,n 3 ),G(ni,n 2 ) = Tpf(ni,n 2 ,n 3 + 1) we have bilinear forms (|50|),(pl|) 
of discrete potential mKdV equation. Thus through the variable transformation (J^j, ((H) the algebraic 
solution of discrete potential mKdV equation (^lj) is given by 



(78) 



P N (m,n-.n 3 ) P N+1 (m,n;n 3 ) 
P N - 2 (m,n;n 3 ) P N ^i(m,n;n 3 ) 

P-N+2(m,n]n 3 ) P- N+ i(m,n;n 3 ) 



P2N-i(m,n;n 3 ) 
P2N -3(111, n;n 3 ) 

Pi(m,n;n 3 ) 



(79) 
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P N (m, n;n 3 + 1) P N+ i(m, n; n 3 + 1) 
fjv-2(m, n;n 3 + 1) P N -i(m, n; n 3 + 1) 



P2N-i(m,n; n 3 + 1) 
-P 2 jv- 3 (m,n; n 3 + 1) 



P-AT +2 (m, n;n 3 + 1) P_Ar +1 (m, n; n 3 + 1) ••• Pi(m,n;n 3 + 1) 

oo 

]T P fe (m, n; n 3 )A fc = (1 - A fll ) m +"(1 - Aa 2 )""(l - A 2 c 2 )-^. 



Concretely, 

P (m, n;n 3 ) 
Pi(m,n;n 3 ) 
P 2 {m,n;n 3 ) 



1, 

aim + (ai — a2)n 
1 



(80) 



(81) 



^a 2 m(m - 1) + 01(01 - a 2 )mn + ^(ai - a 2 ) 2 n 2 - ^(ai + a 2 )(ai - a 2 )n - ^c 2 n 3 , 



where n 3 is an arbitrary constant. 



9 Conclusion 

In this article, we have derived dependent variable transformations for discrete soliton equations to reduce 
them to bilinear forms by using singularity confinement. We have also constructed the A^-soliton solution 
of discrete mKdV equation. We expect that this method based on SC is effective for study on other 
discrete nonlinear equations. 

One of the authors (K.K) was supported by the Grant-in-aid for Encouragement of Young Scientist, 
The Ministry of Education, Science, Sports and Culture, No. 08750090. 

References 

[1] R. Hirota: J. Phys. Soc. Jpn. 43 (1977) 1424. 
[2] R. Hirota: J. Phys. Soc. Jpn. 43 (1977) 2074. 



13 



[3] A. Ramani, B. Grammaticos and J. Satsuma: J. Phys. A: Math. Gen. 28 (1995) 4655. 

[4] B. Grammaticos, A. Ramani and V. G. Papageorgiou: Phys. Rev. Lett. 67 (1991) 1825. 

[5] R. Hirota and S. Tsujimoto: J. Phys. Soc. Jpn. 64 (1995) 3125. 

[6] R. Hirota and S. Tsujimoto: RIMS Kokyuroku 868 (1994) 31. 

[7] S. Tsujimoto and R.Hirota: RIMS Kokyuroku 933 (1995) 105. 

[8] F. W. Nijhoff and V. G. Papageorgiou: Phys. Lett. A 153 (1991) 337. 

[9] T. Miwa: Proc. Jpn. Acad. A 58 (1982) 9. 

[10] A. Ramani, B. Grammaticos and J. Satsuma: Phys. Lett. A 169 (1992) 323. 

[11] Y. Ohta, R. Hirota, S.Tsujimoto and T. Imai: J. Phys. Soc. Jpn., 62 (1993) 1872. 



14 



